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ABSTRACT 

The Type Ib/c supernova SN 2001em was observed to have strong radio, 
X-ray, and Ha emission at an age of ~ 2.5 yr. Although the radio and X- 
ray emission have been attributed to an off-axis gamma-ray burst, we model 
the emission as the interaction of normal SN Ib/c ejecta with a dense, massive 
(~ 3 Mq) circumstellar shell at a distance ~ 7 x 10^^ cm. We investigate two 
models, in which the circumstellar shell has or has not been overtaken by the 
forward shock at the time of the X-ray observation. The circumstellar shell was 
presumably formed by vigorous mass loss with a rate ~ (2 — 10) x 10^'^ Mq yr~^ 
at ~ (1 — 2) X 10^ yr prior to the supernova explosion. The hydrogen envelope 
was completely lost, and subsequently was swept up and accelerated by the fast 
wind of the presupernova star up to a velocity of 30 — 50 km s~^. Although 
interaction with the shell can explain most of the late emission properties of SN 
2001em, we need to invoke clumping of the gas to explain the low absorption at 
X-ray and radio wavelengths. 

Subject headings: stars: mass-loss — supernovae: general — supernovae: indi- 
vidual (SN 2001em) 



1. INTRODUCTION 



The supernova SN 2001em was discovered on 2001 September 20 (Papenkova & Li 2001) 
the galaxy UGC 11794 {z = 0.01935). With an apparent magnitude of about 18.6 and 
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absolute magnitude of M — 16 at the time of discovery (for D = 83 Mpc), the supernova 

was present in an unfiltered image on September 15 but not on September 5 at a level of 
19.5 mag (Papenkova & Li 2001). The spectrum on 2001 October 20 was of Type lb or Ic a 
month after maximum brightness (Filippenko & Chornok 2001). There is little doubt that 
SN 2001em was discovered early, probably before maximum light. We adopt 2001 September 
10, i.e. JD=2452163, as the time of explosion. The 5 day uncertainty in the age is of no 
consequence for the interpretation of the late observations considered here. 

Two years after the explosion, on 2003 October 17, SN 2001em was detected at radio 
wavelengths with the VLA at 3.6 cm. Between then and 2004 January 30, the flux at 3.6 cm 
increased by a factor of ~ 1.3 (Stockdale et al. 2004), and by a factor of 1.56 to 2004 July 1 
(Stockdale et al. 2005). The spectrum in the range 2 — 6.2 cm was nonthermal, Fi, oc 
(Stockdale et al. 2004). The high radio luminosity ~ 2 x 10^^ erg s"^ Hz~^ at 6 cm (Stockdale 
et al. 2004) is unprecedented for a SN Ib/c at this age. Moreover, X-ray observations on 
2004 April 4 (day 937) with Chandra revealed X-ray emission in the 0.5 — 8 keV band with 
a luminosity of ~ 10^^ erg s~^ (Pooley & Lewin 2004), again unprecedented for a SN Ib/c 
at this age. 

The circumstellar (CS) medium in the immediate vicinity of SNe Ib/c is presumably 

shaped by a Wolf-Rayet (WR) type wind. The interaction with this wind is normally ex- 
pected to produce early radio emission that decays with a fairly steep power law in time (e.g., 
SN 1983N, Weiler et al. 1986). Although the emission would normally be undetectable at an 
age of 2 years, Stockdale et al. (2004) undertook late observations of SNe Ib/c based on the 
suggestion that strong late radio emission could be caused by the interaction of a misaligned 
relativistic jet with the WR wind (e.g., Paczyhski 2001). Granot & Ramirez- Ruiz (2004) 
discussed misahgned jet and CS interaction models for the late emission from SN 2001em 
and concluded that the misaligned jet model was favored. A prediction of their model was 
that the radio source should be resolvable with VLBI (very long baseline interferometry) 
observations, with an angular size > 2 mas. VLBI observations have been undertaken by 
2 groups (Stockdale et al. 2005; Bietenholz & Bartel 2005), both of which failed to resolve 
the source at the predicted size; Bietenholz & Bartel set an upper limit on the major axis 
angular size of 0.59 mas. Although the VLBI observations do not support the misaligned jet 
model, the model might not be ruled out if there is relativistic motion of a compact radio 
source. 

However, detection of a strong Ha emission line with FWHM of 1800 km s~^ on 2004 
May 7 (Soderberg, Gal- Yam, & Kulkarni 2004) is not readily explained in the misaligned jet 
model. This type of emission is observed from supernovae which are undergoing strong CS 
interaction and usually associated with Type Iln supernovae (e.g., Filippenko 1997). The 
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high radio and X-ray luminosities of SN 2001em at an age ~ 10^ days are comparable to 
those of bright Type Iln supernovae (e.g., SN 1986J, SN 1988Z) at a similar age. Given this 
similarity, we propose a model of SN 2001em that accounts for the observed phenomena by 
interaction of a SN Ib/c with a dense CS shell that is some distance from the progenitor and 
was initially part of the H envelope of the progenitor star. 

General aspects of our model arc discussed in § 2 and more details are in § 3. The 
model results are compared to observations of SN 200 lem in § 4 and the formation of the 
circumstcllar shell is treated in § 5. The conclusions, with attention to the evolutionary 
status of the SN 200 lem progenitor, are in § 6. 

2. GENERAL CONSIDERATIONS 

We assume that SN 2001em exploded as an ordinary SN Ib/c, and the late X-ray, radio 
and Ha emission were the outcome of the ejecta interaction with a dense CS hydrogen 
shell lost by the progenitor star. This conjecture implies that the CS environment around 
SN 2001em was shaped by two episodes of mass loss: heavy mass loss in a slow red supergiant 
wind (possibly in a common envelope phase or a superwind phase) and a subsequent rarefied, 
fast WR wind. This WR wind caused the formation of a dense CS shell by the interaction 
of the fast wind with a slow dense wind (e.g., Kahn 1983). 

The WR stage following the loss of the hydrogen envelope was relatively brief, so the 
stellar mass could not have decreased significantly at this stage. We expect, therefore, that 
the SN 2001em presupernova was the He core of a massive star. A related example is 
SN 1993 J, which originated from a "typical" mass range 13 — 16 Mq and got rid of almost 
all the H envelope (Woosley et al. 1994). For SN 2001em, we adopt an ejecta mass of 
M ^ 2.5 Mq and an energy of ~ 1.6 x 10^^ erg, similar to parameters found for SN 1993 J 
(e.g., Utrobin 1994); in § 4.1, we consider variation of these parameters. With the adopted 
parameters, the typical velocity of SN material is fsn = (2E/M)i/2 ^ S X 10^ cm s . 
Assuming that X-rays from SN 2001em at t ^ 937 d detected by Chandra (Pooley & Lewin 
2004) correspond to the stage of strong interaction with the CS shell, the radius of the CS 
shell is then Res ~ Vsnt ~ 6 x 10^^ cm. The typical supernova velocity (or less) must be 
chosen if the supernova energy has substantially thermalized; on the other hand, the rising 
radio fiux implies that the interaction is not highly evolved. The hydrogen envelope was 
presumably lost during the red supergiant stage with a velocity ~ 10^ cm s~^. The CS shell 
was likely accelerated by the WR wind to a velocity Ucs ~ 20 km s~^, so the age may be 
Rju,s - 103 yr (§5). 
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The shock interaction of freely expanding SN ejecta {v — r/t) with a CS shell proceeds 
through the formation of a double shock interface layer with the forward shock accelerating 
the CS gas and the reverse shock decelerating SN cjccta. The swept up hot gas between the 
two shocks is responsible for X-rays, while the accelerated electrons and amplified magnetic 
field in the interaction zone bring about the synchrotron radio emission (Chevalier 1982; 
Chevalier & Liang 1989). The Ha emission line with a full width at half maximum (FWHM) 
of ~ 1800 km detected on day 970 (Soderberg, Gal- Yam, & Kulkarni 2004) is probably 
emitted by the CS gas accelerated in the forward shock, since the SN ejecta are devoid of 
hydrogen. The Ha profile has no apparent extended wings beyond a velocity of 2000 km s"^ 
(Soderberg, Gal- Yam, & Kulkarni 2004), which imphes that the bulk of the hne-emitting 
gas in the forward shock moves with velocities < 2000 km s~^. 

Two possibilities for the origin of the high velocity hydrogen in the forward shock are: 
(i) a shock wave with a velocity ~ 2 x 10^ km s~^ that passed through the smooth dense 
CS shell was radiative, so a cool dense shell formed between the shock wave and the contact 
surface; (ii) the dense CS shell was clumpy, so the clouds were first shocked by slow radiative 
shocks, and then fragments of the shocked clouds were accelerated in the forward shock up 
to ~ 2 X 10^ km s^^, similar to the scenario proposed for the Ha emission in SN 2002ic 
(Chugai, Chevalier, & Lundqvist 2004). 

The temperature of the forward shock with a velocity of ~ 2 x 10'^ km s~^ is about 
5 keV, far below the 80 kcV estimated from Chandra observations on day 940 (Pooley & 
Lewin 2004). The reverse shock can be much hotter. As a result of the sudden collision of the 
rarefied outer SN layers with the dense CS shell, the swept up shell at the SN/CS interface 
is strongly decelerated because of a high CS/SN density contrast Pcs/Psn- The swept up shell 
velocity 

^^s ~ Vsn{psn/ Pcsf^ , (1) 

where v^n is the velocity of SN ejecta at the reverse shock, is low and, consequently, the 
reverse shock velocity Vgn — Vs is high enough to provide a high temperature of the shocked 
SN ejecta. An additional factor that favors a high temperature at the reverse shock is an 
average molecular weight of ejecta that is larger than for normal cosmic abundances. For the 
ionized He composition of SN Ib/c, jl = p/nrrip = 1.33, so the velocity of the reverse shock 
must be only ^ 5500 km s^^ to yield a shock temperature of 80 keV. With the velocity of the 
swept up shell ^ 2000 km s^^, the boundary ejecta velocity should be ^ 8000 km s^^ at the 
time of X-ray observation. The actual picture may be more complicated, because the X-ray 
emission is a combination of the radiation from both shocks with different temperatures. 

We now estimate the mass of the shocked SN ejecta assuming that the reverse shock 
is the dominant component. With the standard bremsstrahlung cooling function A = 1.6 x 



- 5 - 



-|^Q-27j'i/2^2^ where Z — 2 is the ion charge for pure He, we find that the observed X-rays 
with T = 9 X 10^ K (i.e. 80 keV) and luminosity ~ lO^-*^ erg s~^ require an emission 
measure EM 5 x 10^^ cm^^. To estimate the mass of the shocked ejecta, the volume 
of the reverse postshock layer between the reverse shock and the contact surface should be 
determined. Let SN ejecta with a power low density distribution p = Po{v/vq)^^ for v > vq 
and p = pq aX V < vq collide with a dense narrow CS shell. The ratio of the radius of the 
reverse shock to the radius of the contact surface ^ = r^s/ Rc in that case is (Chevalier & 
Liang 1989) 



Substituting k ^ 9 gives ^ ^ 0.91. For i?c ~ 6 x 10^^ cm, ^ = 0.91, and EM ^ 5 x 10^^ 
cm~^, the mass of the shocked ejecta is ~ 0.8 Mq. Momentum conservation imphes that the 
CS shell mass required to decelerate Mi ~ 1 Mq of the SN ejecta with the typical velocity 
of Vsn ~ 8000 km s"^ down to ~ 2000 km s"^ is M2 ^ Mi{v^,Jv^ - 1) 3 Mq. Given 
the age of the CS shell of ~ 10^ yr, the last episode of the hydrogen envelope removal thus 
occurred with a mass loss rate of ~ 3 x 10^^ Mq yr~^. Below we present a more detailed 
interaction model that will confirm this general picture, although with some modifications 
regarding the interpretation of the X-ray spectrum (§4.1). 

This consideration of the X-ray luminosity shows that the mass of the CS shell must be 
comparable to the ejecta mass, so that the supernova energy is substantially thermalized. 
Although the radio luminosity was rising at an age of 770-1000 days, we expect that the 
luminosity cannot rise much further because the interaction energy cannot keep increasing. 
This is consistent with the fact that the luminosity of SN 2001em is comparable to that of 
the most luminous radio and X-ray supernovae at a comparable age. The situation can be 
compared to that in SN 1987A, which is known to have a dense ring at a radius of 6 x 10^^ 
cm, about an order of magnitude larger than the dense shell inferred here. The radio flux 
from SN 1987A started to rise at an age of 3 yr (Manchester et al. 2002), but at a luminosity 
level orders of magnitude smaller that the luminosity observed in SN 2001em; the initial 
interaction is with the outer, high velocity, low density ejecta, which have only a small part 
of the supernova energy. In addition, the mass in the CS shell is probably considerably 
less than the ejecta mass, so that only a fraction of the supernova energy will have been 
thermalized at the time that the CS shell is swept up. 




(2) 



-6- 



3. THE MODEL 

3.1. Hydro dynamic Interaction 

We consider freely expanding {v = r/t) SN ejecta interacting with a CS environment 
consisting of a massive dense CS shell placed between the WR wind (inner zone) and the 
slow red supergiant wind (outer zone). The SN ejecta mass is taken to be M 2.5 Mq, 
the kinetic energy is £^ ~ 1.6 x 10^^ erg, and the density distribution has a flat inner zone, 
p — Pq V < Vq, and a power law outer layer, p — Pq{v/vq)~^ for v > Vq with k — 9. 

The WR wind is set by the wind density parameter wi — Mi/ui — 6.3 x 10^^ g cm~^ 
which corresponds to the choice Mi = 10"^ Mq yr~^ and the wind velocity of ui — 1000 km 
s""*^. This wind has little effect on the late time interaction because of its low density, and 
we neglect the fact that the WR wind has to pass through a shock wave inside the CS shell. 
The assumed width of the dense CS shell is R2 — Ri = AR ^ O.IR, where R is the average 
radius of the CS shell. The density in the CS shell is characterized by the density parameter 
W2, which is defined by the shell radius and the shell mass with the assumed AR/R ratio. In 
§ 5.2, we show that the assumed shell thickness and density are consistent with expectations. 
The outer wind at r > i?2 is presumably the slow dense wind of a red supergiant and we take 
^3 = 10^^ g cm~^, an intermediate value between the wind densities expected for a SN IIP 
and a bright SN IIL. Wc explored other values and found that, provided < 10^^ g cm~^, 
the results are not sensitive to this parameter. The two versions of the density distribution 
of the CS matter (CSM) that will be used below are shown in Fig. 1. The models differ in 
the radius of the CS shell. In model A, the interaction with the CS shell is not yet finished 
by the time of the X-ray observation (t = 937 d), while model B, with a smaller shell radius, 
represents a situation where the CS shell was swept up prior to the X-ray observation. 

The hydrodynamic evolution of the SN interaction is treated in the thin shell approx- 
imation (Chevalier 1982), i.e. the layer between forward and reverse shocks is replaced by 
a thin layer. This description yields the radius of the thin shell Rg, the shell velocity Vg, 
and the boundary velocity of the unshocked SN ejecta Vsn at an age t. The velocities of 
the forward shock (f« Vg) and reverse shock {vsn — Vs) provide us with kinetic luminosities 
for both shocks. The kinetic luminosities can be converted into X-ray luminosities if they 
are multiplied by a factor te/ {te + tc) which is determined by the ratio of the cooling time 
tc and the time of accumulation of the shock internal energy tc = {d\nE/dt)-^. Usually 
te = t is assumed (e.g.. Chevalier & Fransson 1994; Chugai 1992) and this is appropriate for 
interaction with a smooth wind. In the case of a collision with a narrow CS shell one expects 
te < t, so it is more appropriate to use the directly calculated value of t^- However, this 
approach based upon the instant kinetic luminosity is not valid at a very late epoch, when 
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the forward shock has overrun the dense shell {Rs > R2) and propagates into the rarefied 
wind. The radiation from the large mass of shocked CS gas left behind by the previous 
interaction should dominate. We treat the thermal history of the outer hot gas by solv- 
ing the time- dependent energy balance with the internal energy generation in the forward 
shock, adiabatic cooling (dominant term), and radiative cooling. To obtain the density of 
the outer hot gas, we assume that the shell expands homologously and AR/R = 0.1. This 
approximation is supported by the fact that the spherical expansion of a gas cloud proceeds 
in a self-similar way with a boundary velocity comparable to the initial thermal velocity 
(Zel'dovich & Raizer 1967), i.e. the shell speed in our case. 

The model simulations show that model A assuming a smooth CS density predicts 
a larger absorption of X-rays than indicated by the observations. Therefore, we assumed 
a clumpy structure of the CS shell to reduce the absorption. In model B the X-rays do 
not constrain the structure of the CS shell, because it has been already shocked at the 
time of observation. We approximate the interaction with a clumpy CS shell by including 
the dumpiness only for the absorption computation, while the interaction dynamics and 
the unabsorbcd X-ray luminosity of both shocks are described assuming a smooth CS shell 
with an average density. This approximation, although crude at first glance, is justified if 
CS clouds entering the forward shock fully deposit their kinetic energy, which means that 
clouds are completely crushed, fragmented and mixed within the forward shock. This picture 
amounts to the assumption that, for CS clouds of radius a, the cloud crushing time tec — a/vc 
(Klein, McKee, & Colella 1994) is significantly less than the time required for a CS cloud to 
cross the forward shock region with width h, i.e., ^cc ^ {h/R)R/vs. We return to this issue 
below (§4.1). 

The question arises of whether model A with a smooth density can adequately represent 
the X-ray spectrum. The spectrum of the forward shock in the cloudy CS gas is a combination 
of radiation from the hot intercloud gas and cooler cloud shocks, which cannot be treated 
properly in the smooth approximation. Fortunately, the X-ray emission of cloud shocks does 
not affect markedly the amount of internal energy generated by the forward shock in the 
intercloud gas. The point is that the total kinetic luminosity of slow cloud shocks is a factor 
of Xi -C 1 (where Xi 1 is cloud-to-intercloud density ratio in the forward shock) smaller 
than the kinetic luminosity of fast bow shocks. Therefore, most of the internal energy in 
the forward shock region resides in the hot intercloud gas. The approximation of a smooth 
forward shock region for the clumpy model A may be satisfactory for the X-ray spectrum in 
the range E > 1 keV. Yet a sizable amount of the kinetic luminosity of the forward shock 
interacting with a clumpy CS shell is expected to be emitted by cloud shocks in the soft 
band {E < 0.5 keV). Most of this radiation is presumably absorbed and re-emitted in the 
optical/UV band outside of the X-ray band of interest. 
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3.2. X-ray and Radio Absorption Model 



The observed X-ray spectrum is affected by four cool components, including (1) galactic 
absorption and intrinsic absorption in (2) unshocked SN ejecta; (3) unshocked CS gas; and 
(4) shocked cool Ha-emitting gas. The reddening towards SN 2001em is E{B — V) = 0.1 
(Schlegel, Finkbeiner, & Davis 1998), which corresponds to a column density A'h = 6 x 
10^° cm"^ (Spitzer 1978). We take this value as the total column density in the first cool 
component, neglecting interstellar absorption in the host galaxy. The SN ejecta can absorb 
X-rays, despite the fact that ejecta gas can be heated up to ~ 6 x 10^ K. The smooth outer 
wind in the region r > R2 irradiated by X-rays becomes hot (T > 10® K) for U73 < 10^® 
g cm~^ and is not a significant absorber of radio. However the outer wind absorbs X-rays, 
although carbon and oxygen may be ionized up to He- and H-like ions. 

As noted above, the CS shell in model A has to be clumpy to be consistent with the 
low X-ray absorption. The optical depth of the clumpy CS shell is computed as 

f-Rb 



where i?b is the outer radius of the CS region taken to be 1.5 x 10^'' cm, Ny = f /{Ana^ /3) 
is the number density of clouds, / = (1 — {)/% is the cloud filling factor, ^ is the mass 
fraction of the intercloud phase (we take ^ = 0.1) and x is the ratio of cloud density to 
the local average density, Tc — {4:/3)ak{E)pc is the cloud average optical depth, and k{E) 
is the absorption coefficient for X-rays with energy E. The occultation optical depth of 
the CS shell along the radius, Tqc = TTa^NvA.R, should be of order unity or less to provide 
reasonable transparency. With fixed and x, the only remaining parameter that determines 
Toe is the cloud radius a. The X-ray absorption by the unshocked SN ejecta is determined 
by integrating the radiation transfer equation along rays with different impact parameters 
for the assumed density distribution. We found that the absorption of radio emission by 
SN material has a small effect on the emergent spectrum, so we treat it approximately in 
terms of the average values of the electron number density and the temperature of the ejecta 
irradiated by X-rays. 

In model B, only the outer wind and fourth component (Ha-emitting gas) can provide 
intrinsic X-ray absorption. To treat the latter we assume that 10 — 30% of the shocked CS 
gas resides in the cool Ha-emitting gas; the result is insensitive to the value in this range. 
This gas is probably distributed in the form of heterogeneous structures (knots, filaments, 
sheets) imbedded in the hot gas of the forward shock. We assume that hot and cool gas are 
homogeneously mixed. In this case the absorption by the fourth component is characterized 
by the optical depth T4,av for the average column density, and by the occultation optical 
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depth (the average number of clouds along the radius) r^^oc- The optical depth for the X-ray 
radiation transmitted through the forward shock layer derived from Eq. (3) is then 

T4 = 'r4,oc[l - exp (-r4,av/r4,oc)]- (4) 

For photons emitted in a layer along the normal to the surface, the intensity of escaping 
radiation is / = /o[l — exp (— r4)]/r4, where the /q is the unabsorbed intensity. We adopt 
the same expression for the flux of radiation escaping from layer F = Fo[l — exp (— r4)]/r4. 
However, for the radiation escaping through the inner surface we additionally have to take 
into account the absorption in the SN ejecta (e^^^) and absorption in the Ha-emitting gas 
at the opposite side of the forward shock (e""^*). For the radio emission, we assume also that 
the Ha-emitting material is homogeneously mixed with the radio-emitting shell. This is an 
approximation to a complicated situation in which the forward and reverse shock regions 
both contribute to the acceleration of relativistic electrons and magnetic field amplification, 
so, generally, the radio-emitting shell and the Ha-emitting material may not overlap exactly. 

The free-free absorption of the radio emission by the CS shell depends on the ionization 
fraction of hydrogen (x) and the electron temperature of the cool components. These values 
are calculated taking into account that the absorbed energy of X-rays is shared between 
Coulomb heating, excitation and ionization. We assume that the latter two processes have 
equal branching ratios. For the Coulomb heating we adopt a branching ratio , a rea- 
sonable approximation if the ionization fraction x > 10~^ (Kozma & Fransson 1992). The 
cooling term in the energy balance is calculated using a standard cooling function for solar 
composition (Sutherland & Dopita 1993). If the thermal balance cannot maintain the cool 
phase (T < 10^ K) we adopt full ionization and consider this gas as transparent. The typical 
value of the ionization fraction of the clouds at about 900 day is x ~ 0.3 and temperature 
Te ^ 13000 - 14000 K, while the intercloud gas is hot {T^ > 10^ K) and thus does not absorb 
radio emission. For the Ha-emitting gas the typical ionization is a; ~ 0.01 and temperature 
Tc ^ 10^ K. In the SN ejecta the helium is fully ionized and the electron temperature is 
Te «i 6 X 10^ K. 

The unabsorbed X-ray spectrum is modelled by thermal bremsstrahlung with a non- 
relativistic Gaunt factor (Itoh et al. 2000). The photoionization cross sections for the ab- 
sorption by electrons of K and L shells with 2s and 21 subshells are taken from Verner 
et al. (1993). Solar abundances are assumed for the CSM, while the SN Ib/c ejecta are 
approximated by a mixture of 65% He, 33% O, and 2% Fe by mass . 
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4. RESULTS 

4.1. CS Interaction Models and X-ray Emission 

We now discuss models that fit the basic observations of SN 2001em. The parameters 
of the CS shell for models A and B, namely, the inner and outer radius of the CS shell {Ri 
and i?2), the density parameter of the CS shell (u»2), and the total mass of the CS shell are 
given in Table 1 (see also Fig. 1). The total mass of the CS shell is 2.3 - 3 Mq (Table 1), 
similar to the preliminary estimate in § 2. The evolution of the major output parameters is 
presented in Fig. 2. The plot shows the radius of the thin shell (contact surface), the velocity 
of the thin shell and of the boundary of the SN cjccta, the unabsorbcd X-ray luminosity, and 
the temperature of forward and reverse shocks. Here we assume a smooth CS shell so both 
shocks are adiabatic; however, the forward shock between days 400 and 500 is very close to 
the radiative regime with tc/U ~ 1.5 — 2. 

After about day 200 a shell formed during the interaction of the SN with the WR 
wind coUides with the inner boundary of the dense CS shell. There is subsequent rapid 
deceleration from ~ 25000 km s~^ down to 1200 — 1300 km s~^, followed by a period of 
steady acceleration (e.g., Chevaher & Liang 1989; Dwarkadas 2005). The collision results in 
a rise of the X-ray luminosities of both shocks toward their maximum values, ~ 10^^ erg s~^. 
The maxima occur at the phase when the CS shell has been swept up, i.e. at ^ 1000 d and 
^ 900 d for models A and B respectively. The luminosity of the shocked CS gas remains 
high after the CS shell has been overtaken; the temperature, however, rapidly decreases for 
t > 1000 d due to adiabatic cooling. The contribution of the forward shock to the internal 
energy of the shocked CS gas at a late epoch is negligible for t < 2000 d. At the time of the 
Ha observation {t — 970 d), the thin shell velocity is 2200 km s"-*^ in model A and 3200 km 
s~^ in model B; both are consistent with the observed Ha width. 

The calculated X-ray spectra for both models (Fig. 3) reproduce the data of Pooley 
& Lewin (2004) quite satisfactorily. Note that model A has a clumpy CS shell here; a 
smooth CS shell produces too strong absorption (Fig. 3). The fit quahty for both models is 
comparable to that of the isothermal hot spectrum with kT — 80 keV (Fig. 3, inset) and with 
external absorption corresponding to A^h = 1-6 x 10^^ cm~^, the value reported by Pooley & 
Lewin (2004). The unabsorbed spectrum in our models is a mixture of the hot (kT ~ 100 
keV) radiation of the reverse shock and "cool" (kT ~ 5 — 6 keV) radiation of the forward 
shock with comparable luminosities (Fig. 2). In fact, the cooler component dominates in 
the range < 5 kcV. We thus come to a picture in which the data are reproduced by a 
relatively cool spectrum subjected to the intrinsic absorption by the clumpy cool material, 
so the emergent spectrum succesfuUy mimics the observed hot isothermal spectrum (~ 80 
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keV) with a smooth external absorber. Although the general picture presented earlier (§ 2) 
is based on the view that the X-rays arc represented by a hot isothermal spectrum, the 
mass estimate is not appreciably changed, because the luminosity of the hot component is 
comparable to the total luminosity. 

The Ha-emitting gas is characterized by an occultation optical depth Toc,4 = 0.5 in 
model A and roc.4 = 1-5 in model B. The cloudy CS shell in model A has Toc,2 = 1- We 
adopt a cloud density contrast (cloud-to-average) x = IOj which leads to a cloud radius 
a = 4.7 X 10^"^ cm. The result is not sensitive to variation of x in the range 3 — 30. With 
these parameters we can check the condition for the approximation of a smooth density in 
model A, i.e. a <C hvc/vg = ao (§ 3.1). From momentum conservation, the ratio of the cloud 
shock velocity to the forward shock velocity Vc/vs ~ ~ 0.2. Assuming the forward 

shock width h ~ O.IR Ri 7 x 10^^ cm, this gives the upper limit ao = 1.4 x 10^^ cm which 
is three times larger than the cloud radius a — 4.7 x 10^^ cm in model A. The condition for 
the smooth density approximation is thus barely met. 

To consider variations in the supernova parameters, we explored "low energy," M — 
2.5 Mq, E = 10^^ erg, and "high mass," M = 4 M©, £^ = 1.6 x 10^^ erg, versions of model 
A. We found that the X-ray spectrum is well reproduced in both cases, if the CS shell radius 
is smaller by a factor 0.8 — 0.84. The thin shell velocity is 1800 km s~^ in the low energy 
case and 2000 km s~^ in the high mass case, i.e. still consistent with the Ha line profile. 
The uncertainty in the SN ejecta parameters thus results in minor variations in the model 
parameters. 

4.2. Radio Emission 

The collision of the SN ejecta with the CS environment is expected to produce syn- 
chrotron radio emission as a result of the rclativistic particle acceleration and magnetic field 
amplification in the shock wave region (Chevalier 1982). As noted above, the situation with 
the radio evolution of SN 2001em is similar to SN 1987A, in which the radio fiux showed 
steady growth after the forward shock began to interact with the dense H II region in the red 
supergiant wind (Chevalier & Dwarkadas 1995; Manchester et al. 2002). Here we concentrate 
on the absorption effects of both free-free absorption and synchrotron self-absorption, and 
consider some implications of the radio evolution. Our analysis is based on the assumptions 
that rclativistic electrons with a power law spectrum, dN/dE = KE~^ (the spectral index 
is a = (p — l)/2 for oc z/~"), and magnetic field B are distributed homogeneously in a 
spherical shell with width Ai? = 0.1i?s and radius 7?s; the radius is provided by the inter- 
action dynamics. The radio-emitting shell is assumed to be homogeneously mixed with the 
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cool clumpy Ha-emitting gas. We assume a minimum energy of the relativistic electrons of 
1 MeV, energy equipartition between magnetic field and relativistic particles, and equiparti- 
tion between relativistic electrons and protons. Given the interaction model, the remaining 
free parameter is the ratio {() of the energy density of magnetic field plus relativistic particles 
to the kinetic energy density of the forward shock. 

The effects of SSA (synchrotron self-absorption) and free-free absorption on the radio 
spectrum in model A and model B are shown in Fig. 4 along with data on 2004 January 31 
(Stockdale et al. 2004). The optimal parameters are a = 0.55, C = 0.024 in model A and 
a = 0.5, ( = 0.017 in model B. We found that the variation of the SN cjccta parameters in 
the low energy and high mass models (§ 4.1) results in a 10% variation of the parameter ( 
with practically the same fit to the spectrum on day 872. The derived values of ( imply a 
magnetic field B 0.2 G in both models at the time under consideration. For this magnetic 
field, synchrotron losses become important (i.e. tgyn ~ 1 yr) at frequencies u > 2 x 10^^ 
Hz. In the observed frequency range {u < 2 x 10^° Hz), the spectrum is not affected by 
synchrotron cooling. 

The spectrum without free-free absorption demonstrates that SSA is significant for 
i/ < 4 X 10^ Hz. The importance of SSA could be deduced from equation (13) of Chevaher 
(1998), which shows that the size of the radio emitting region at which the shell turns 
optically thin is comparable to the size of the shell found here. An observation of SN 
2001cm on 11 March 2005 at 1.6 GHz showed a considerably lower flux than would have 
been expected from the higher frequency flux evolution (Paragi et al. 2005), which can be 
interpreted as a low frequency turnover or a sudden flux decline. Our results are consistent 
with a low frequency turnover. The spectral inversion at u > 2 x 10^° Hz is related to the 
transparency of absorbing clumps at high frequencies. 

The mechanisms for particle acceleration and magnetic held ampliflcation, especially in 
the case with a cloudy structure of the GS shell, are uncertain. The evolution of the radio 
fiux, therefore, cannot be predicted with confidence. Assuming that the energy of relativistic 
electrons and magnetic field is a constant fraction of the kinetic energy (i.e., ( — const), we 
find that in model A the radio fiux at 3.6 cm on day 767 is ~ 3 times lower compared to 
the observed flux, while on day 1025 the model flux is ~ 2 larger. This implies that the 
parameter ( should drop by factor of three between days 767 and 1025. In model B assuming 
a smooth GS shell, the flux on day 767 is 30 times lower owing to strong absorption in the GS 
shell. On the other hand, assuming a clumpy GS shell with the same dumpiness parameters 
as in model A, we are able to reproduce in model B the flux evolution between day 767 
and 1024 with ( Ki const. This analysis indicates that model B with a clumpy GS shell is 
somewhat preferred from the point of view of the interpretation of X-ray and radio data. 
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4.3. HcK-emitting Gas 



As remarked above (§ 2), the Ha emission line on day 970 (Soderberg, Gal- Yam, & 
Kulkarni 2004) indicates the presence of cool shocked CS gas behind the forward shock 
front. We stress that in model B with the overtaken CS shell, the forward shock is adiabatic 
for a smooth density distribution and the assumed CS width AR/R — 0.1. The shock 
becomes radiative if the density is at least factor of three larger, i.e. n > 2 x 10^ cm~^. This 
requires either a factor of three narrower CS shell, i.e. AR/R < 0.03, or, alternatively, the 
CS shell could be clumpy. The latter option is preferred because the radio data, as noted 
above, also suggest a clumpy structure for the CS shell in model B. To be consistent with the 
Ha observation, the cool gas in the forward shock of model B must survive at least during 
~ 100 days after it has been shocked. This is plausible, although the survival of cool gas in 
the hot environment is a complicated issue, because it is related to the poorly defined mass 
exchange between cool and hot phases, especially in the presence of magnetic field. 

The dense shocked cool hydrogen re-emits the absorbed X-rays, producing Ha. The Ha 
luminosity can be written as a fraction of the X-ray luminosity of the forward shock absorbed 
by the fast cool hydrogen multiplied by the efficiency rj of the Ha emission. The fraction 
of absorbed X-ray radiation with the bremsstrahlung spectrum Fe oc E~'^-^ exp {—E /kT) 
(Cox 2000) is ~ {E-i/kT)^-^, where Ei is the energy at which the optical depth of the 
cool hydrogen shell with a mass Mcooi is unity. An energy dependence of the absorption 
coefficient /Cx(-E') oc E~^/^ implies Ei cx M^/^j. The expected Ha luminosity on day 970 with 
the temperature of X-ray radiation of the forward shock kT = 5 keV (Fig. 2) is then 



where I/x,4i is the X-ray luminosity of the forward shock in units of 10^"^ erg s~^. The 
efficiency of Ha emission according to photoionization models for relevant parameters, i.e. 
the pressure equilibrium density of cool (~ 10^ K) gas of ~ 3 x 10^^ cm column density 
of ~ 10^^ cm and X-ray flux of ~ 10^ erg s~^ cm~^ is ~ 0.1 (Collin-Souffrin & Dumont 
1989, model 55). The swept up CS mass in the forward shock on day 970 is 2 Mq in model 
A. Assuming one-third of it resides in the cool phase (0.7 Mq) the expected luminosity of Ha 
(Eq. [5]) is then ~ 4.5 x 10^^ erg at about day 1000 for Lx,4i = 1. A comparable amount 
of a cool gas in model B would produce a similar Ha luminosity. These estimates ignore 
the dumpiness of the Ha-emitting gas that might slightly reduce the luminosity. On the 
other hand, we also ignored the effect of thermal coductivity that might somewhat increase 
the Ha luminosity. Unfortunately, the observed Ha flux is not yet available; its value would 
provide a useful constraint on models. 




(5) 



An imphcation of the SN/CSM interaction model is the possible existence of a narrow 
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HcK line from the ionized pre-shock CS gas. For model A the luminosity of the narrow Ha is 
~ 10^^ erg s"-*^ on day 1000, i.e. ~ 2% of the broad Ha line. In model B with the overtaken 
CS shell, a narrow Ha line is not expected. A search for narrow Ha emission could provide 
a test for the presence of dense unshocked CS gas. 

5. FORMATION OF THE CS SHELL 

5.1. Age of the CS Shell 

The large mass and high density of the CS shell suggest the following formation scenario: 
loss of the hydrogen envelope ~ 10^ yr before the supernova, possibly during a common 
envelope phase or the superwind phase of a single star (Heger et al. 1997), and the subsequent 
sweeping up of this matter by the fast wind of the hot He star. We thus envisage three major 
phases of the CS shell formation: (I) mass- loss of the hydrogen envelope {0 < t < ti); (II) 
sweeping up the shell by the WR wind {tj < t < tu); and (HI) acceleration of the swept 
up shell by the WR wind {tu < t < tm). Garcia-Segura et al. (1996) have simulated this 
sequence of events, although for different stellar parameters. To estimate the age at the 
moment of the SN explosion {tm), we use the average parameters of our models: final CS 
shell radius rm = 6.5 x 10^^ cm and the mass of the CS shell Mcs = 2.7 Mq. The velocity 
of the mass-loss at stage I is assumed equal to the typical red supergiant escape velocity 
Ursg — 10 km s~^, while the mass- loss rate at this stage Mrsg is a free parameter to be 
determined. For the WR wind, the assumed parameters are the same as stated above, i.e. 
u^J^^ = 1000 km s~^ and Mwr = 10~^ Mq yr'^. 

The mass-loss rate at the first stage can be constrained using the following arguments. 
The minimum value of Mrsg is determined by the condition that at least the two first dy- 
namical phases should pass with the correct final radius of the CS shell (r f« 6.5 x 10^^ cm). 
The second phase is needed to build the density of the CS shell up to a value consistent 
with our models. An upper limit of M^sg ~ 0.01 Mq yr~^ is obtained assuming that the 
average density parameter W2 ~ 6.5 x 10^^ g cm~^ in our models (Table 1) corresponds to 
the unperturbed CS shell ejected by the presupernova. 

The duration of the first phase is ti — Mcs /Mrsg- To estimate the time length of phase 
II, we consider the interaction of two steady winds. The interaction proceeds in the pressure- 
dominated regime, i.e. the WR wind termination shock has a relatively small radius. In that 
case the swept-up shell expands with the constant velocity Atirsg determined by the cubic 
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equation (Kahn 1983) 



1 MwR / "UWR \ 



2 



A(i-A)^ = -^^P^ . (6) 



In phase III, a thin shell of mass M is driven by the pressure of the shocked WR wind 
with the kinetic luminosity I/wr a-nd the internal energy of the shocked wind (bubble) E. A 
similar problem for the shell pushed by a pulsar wind has been solved earlier (Reynolds & 
Chevalier 1984). The difference is in the adiabatic index: here we have 7 = 5/3 instead of 
4/3. The equation of motion of the shell is then 

..d^r 2E 



while the energy equation is 



dE _ 2E dr 

^-i^WR-— ^. (8) 



The solution of these equations is the self-similar evolution of the CS shell radius 

/ 2Lwr\ ^^"^ ^3/2 



1/2 



To apply this solution to the phase III of the CS shell expansion, we identify the period 
til < t < till with the time taken by the CS shell to expand in the self-similar regime from 
the radius rn at t — tu to the final radius rni- 

For the parameters adopted in model B and the maximum mass-loss rate during stage 
I (Mrsg « 0.01 Mq yr"^), the phases I, II, and III end at ti = 260 yr, tu = 480 yr, and 
tin ~ 830 yr respectively. The CS shell velocity at the moment tm is vm — 55 km s~^. The 
existence of the acceleration phase tu < t < tui provides in this case a natural mechanism for 
the fragmentation of the CS shell by the growth of the Rayleigh-Taylor instabihty (Garcia- 
Segura et al. 1996). The minimum mass- loss rate at stage I, determined from the requirement 
that stage II should end just prior to the SN explosion, is Afisg ~ 0.002 Mq yr~^. In this 
case ti — 1400 yr, tu — 2000 yr, and the shell velocity is 33 km s~^. 

We thus conclude that a major mass-loss episode with the rate (2 — 10) x 10~^ Mq yr~^ 
took place 1000-2000 yr before the SN explosion. The WR wind accelerated this shell up to 
velocity of 30-50 km s~^. 



5.2. Density of the CS Shell 

The minimum density in the CS shell predicted by our models is n > 2 x 10^ cm~^ 
(§4.3). This value should reflect the pressure and thermal balance at the latest phase of the 
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interaction of the WR wind with the CS shell. The dynamical pressure created by the WR 
wind is 

For the parameters adopted above (Mwr = 10~^ Mq yr~^, fwR = 1000 km s~^, R = 
6.5 X 10^^ cm) this expression leads to pwR ~ 10^^ dyn cm^^. Because pwR is actually 
determined by the pressure at the WR wind termination shock, the pressure could be a 
factor ~ 30 higher than this estimate if the wind bubble is contained within the shell. 
However, if the bubble can break out of the shell because of instabihties, the pressure is 
reduced. We regard this pressure estimate as a lower limit. To derive the density from 
pressure equilibrium one needs to know the temperature of the compressed CS gas. First, 
we check that the CS shell was mostly neutral. 

The ionizing radiation of the presupernova forms a low density HII zone. The density 
of the HII zone in pressure equilibrium is n2 — PwR/kT2 ~ 7 x 10^ cm~^ assuming T2 = 10^ 
K. The column density of the HII zone is 

N2 = ^^ , (11) 

where Luv is the luminosity of the ionizing radiation, a is the hydrogen recombination 
coefficient, and Ry is the hydrogen ionization potential. The luminosity of the progenitor, 
assuming a hehum core of a 15 Mq star, is L = 2 x 10^^ erg s~^ (Meynet et al. 1994) . Adopting 
Luv = 0.5L ^ 1 X 10^8 erg s-\ T2 = 10^ K, a = 2 x lO'^^ s-\ and = 6.5 x 10^^ 
cm, we obtain A'"2 ~ 6 x 10^° cm~^. This is substantially lower than the hydrogen column 
density of the CS shell, ~ 6 x 10^^ cm~^ for a CS shell with a mass of 2.7 Mq. This estimate 
implies that the presupernova radiation ionizes only a small fraction of the CS shell and the 
bulk of the CS shell is neutral and can be referred to as a photodissociation region (Tielens 
& HoUenbach 1985). 

The cool neutral CS shell with a column density of A^h ~ 6 x 10^^ cm^^ is characterized 
by a dust extinction Ay ~ 10, assuming interstellar dust properties (Spitzer 1978). In a 
photodissociation zone with dust extinction Ay > 4 the gas temperature follows the dust 
temperature owing to the gas heating by the dust IR radiation (Tielens & HoUenbach 1985). 
With the emission efficiency of the dust at the maximum of the blackbody infrared emission 
Qe — a/X ^ 3.4aTd (here a is the grain radius), the temperature of the dust grains should 
be Td ^ {L/557rR^aay/^. With L = 2 x 10^^ erg s'^ R = 6.5 x 10^^ cm and a = 10"^ cm, 
we obtain 300 K. Gas with this temperature in equilibrium with the pressure pwR has 
a density n ~ 2 x 10^ cm~^. This is consistent with the lower limit for the density of CS 
shell clouds. The suggested scenario of the CS shell formation thus seems realistic. 



-17- 



6. DISCUSSION AND CONCLUSIONS 

Our goal was to propose a model for the strong late time X-ray, radio, and Ha emission 
from the Type Ib/c supernova SN 2001em. We developed a picture in which the SN ejecta 
of a normal SN Ib/c collide with a dense massive CS shell (M^s ~ 3 Mq) at a radius of 
~ 0.7 X 10^^ cm. We found that two scenarios are viable: model A in which the dense CS 
shell not is yet overtaken by the time of X-ray observation (937 d) and model B with the 
CS shell overtaken prior to the X-ray observation. The mass loss episode that led to the 
formation of the CS shell was characterized by a high mass loss rate (2 — 10) x 10~^ Mq yr~^ 
at a stage 1000-2000 yr prior to the SN outburst. This material was then swept up by the 
fast WR wind from the presupernova progenitor. Rayleigh- Taylor instabilities likely led to 
the fragmentation of the CS shell which is manifested in the low X-ray and radio absorption 
despite the high column density of the shell. 

The mechanism for the rapid loss of the 3 Mq hydrogen envelope is not clear, although 
it is probably related to mechanisms for producing SN Ib/c presupernovae in a close binary 
scenario. The bulk of SN Ib/c presupernovae are thought to originate by the stripping of 
the hydrogen envelope at a late stage of binary star evolution as a result of mass transfer 
to a companion and/or the loss of the common envelope (van den Heuvel 1983; Ensman & 
Woosley 1988; Podsiadlowski, Joss, & Hsu 1992; Wellstein & Langer 1999). In the case of 
SN 2001em, we believe that common envelope evolution is the most likely mechanism of the 
formation of the massive CS shell. The 1000 — 2000 yr delay between a strong mass loss 
episode and the supernova explosion indicates that the former might happen at the stage of 
carbon burning in the core of a star with an initial mass in the range 14 — 17 Mq (Hirshi, 
Meynet, & Maeder 2004). We suggest that the remains of the hydrogen envelope (~ 3 Mq) 
were lost during the carbon burning phase due to the formation of a common envelope in a 
binary system. 

The relative fraction of SNe Ib/c passing through the common envelope phase is uncer- 
tain and lies between 10% (Podsiadlowski, Joss, & Hsu 1992) and nearly 100% (Tutukov & 
Yungelson 2002). The massive star acquires a red supergiant structure at the He burning 
stage. Therefore, most SN Ib/c presupernovae pass through the common envelope stage long 
before (~ 10^ yr) the supernova explodes in the rarefied CS environment formed by the WR 
wind. Only a small fraction of SN Ib/c presupernovae lose their H envelope at the C burning 
stage. It was argued that this fraction is possibly about ~ 10~^ (Chugai 1997). We thus 
expect that roughly 0.001 — 0.01 of all SNe Ib/c possess close {R < 10^^ cm) massive CS 
shells and display CS interaction by strong radio, optical and X-ray emission at an age of 
1 — 10 yr. Upper limits on late radio emission from SNe Ib/c (Soderberg, Nakar, & Kulkarni 
2005) seem to be consistent with the estimated fraction of these events. 
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The occurrence of the supernova close to the time that the H envelope was lost may 
have relevance to some Type Iln supernovac with strong circumstellar interaction. (Pransson 
et al. 2002) found that the cjecta in SN 1995N inside the reverse shock wave were heavy 
element rich, implying that this Type Iln supernova had lost essentially all the hydrogen 
envelope in a dense wind at the time of the supernova. It is unlikely, however, that the rate 
of all the Type Iln supernovae can be explained entirely by the common envelope events at 
the carbon burning phase (Chugai 1997). 

In addition to radio detection, an efficient way to find SNc Ib/c with close CS shells 
could be the detection of late time Ha emission with a luminosity of ~ 10^^ — 10^° erg s^^ at 
an age of 1 — 10 yr, as has occurred for some Type Iln supernovae. Another interesting way to 
detect the presence of the massive close CS shell around SN Ib/c is related to the possibihty 
that the shell is dusty. If the radius of the CS shell exceeds the dust evaporation radius 
~ (1 — 3) X 10^^ cm, the observations of the infrared echo in the KLM bands may reveal the 
presence of a CS shell soon after the SN Ib/c explosion, similar to infrared dust cchos observed 
from SN 1979C and SN 1980K (Bode & Evans 1980; Dwek 1983). These observations will 
also serve to distinguish late circumstellar interaction from misaligned gamma ray burst 
events, which must be present in the population of SNe Ib/c at some level. 

We are grateful to Dave Pooley for kindly sending us the observed X-ray fiuxcs of 
SN 2001em. This research was supported in part by Chandra grant TM4-5003X and NSF 
grant AST-0307366. 
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Table 1: Model parameters 



Parameter 


Units 


A 


B 


Ri 


10^^ cm 


6.3 


5.3 


R2 


10^6 cm 


7 


5.9 


W2 


10^^ g cm~^ 


7 


6 




Mo 


3.1 


2.3 
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Fig. 1. — The circumstellar density distribution for model A {thick line) and model B {thin 
line) (see Table 1). Inset: the same distributions in logarithmic coordinates. 
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Fig. 2. — Evolution of the circumstellar interaction of SN 2001em: model A on the left (a) 
and model B on the right (b). The plot shows from top to bottom the radius of the thin shell, 
the expansion velocity of the thin shell [solid) and the SN ejecta velocity at the reverse shock 
(dashed), the unabsorbed X-ray luminosity of the forward shock (solid) and reverse shock 
(dashed), and the temperature of the shocked CS gas (solid) and shocked SN gas (dashed). 
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Fig. 3. — The photon spectrum of escaping X-ray radiation on day 937. Both model A 
{thick line) with a clumpy CS shell and model B {thin line) show satisfactory agreement 
with the observed fluxes (Pooley & Lewin 2004). The models include the absorption by the 
clumpy cool Ha-cmitting gas in the forward shock. The fits are comparable to the fit of an 
isothermal hot gas (T = 80 keV) model with an external cool absorber with A^h = 1-6 x 10^^ 
cm~^ (see inset, Pooley & Lcwin 2004). Model A with a smooth CS shell {dashed line) shows 
unacceptably strong absorption. 
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Fig. 4. — The model radio spectrum compared to the data of Stockdale et ah (2004) (dots). 
The upper panel shows the radio spectrum for model A without free- free absorption {dashed 
line) and with free-free absorption {solid line) produced by a clumpy CS shell and Ha- 
emitting gas in the forward shock. The lower panel is the same as the upper panel but 
for model B and the free-free absorption produced only by Ha-emitting gas in the forward 
shock. 



